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Abstract 

We present a theoretical interpretation of the recently revealed features of tempera- 
ture evolution in the ultracold plasma clouds released from a magneto-optical trap, 
namely: (a) its independence at the sufficiently large times on the initial plasma pa- 
rameters and (b) the asymptotics close to t _1 instead of t~ 2 , expected for a rarefied 
ideal gas. It is shown that both these properties can be well explained by the model 
of virialization of the charged particle velocities in the regime of strong electron-ion 
correlations, while heating due to inelastic processes (e.g. three-body recombina- 
tion) should be of secondary importance. These conclusions are confirmed also by 
the results of ah initio computer simulations. 
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1 Introduction 



A new interesting branch of plasma physics, emerged in the last decade, is 
studying the clouds of rarefied ultracold plasmas produced by the laser capture 
and cooling of gases in the magneto-optical traps (MOT) and their subsequent 
ionization (for general review see, for example, [TP1I5] and references therein). 
These are the classical (non-quantum) gaseous systems with characteristic 
temperatures from a fraction to a few Kelvin, in which the Coulomb's coupling 
parameter Y = e 2 nC 3 /ksT can reach considerable values. For example, the 
measured values of this parameter for the ions I\ are about 2-=-3 [4]; the 
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estimates of T e for the electrons are less accurate and model dependent, but 
they also results in the values comparable to unity. 

The possibility of existence of such metastable plasmas was theoretically pre- 
dicted a quite long time ago (e.g. article [5] and references therein), but they 
were created experimentally only after sufficient development of the technology 
for laser cooling of atoms. Besides, it was proposed to produce similar plasma 
states by the artificial release of the ionized gas clouds from spacecraft [6|7j . 
but the possibility of diagnostics in cosmic space remains too limited by now. 

One of the most interesting recent results in the experimental study of ultra- 
cold plasmas is the behavior of temperature in the clouds released from a trap 
and expanding freely in space. It was unexpectedly found that, firstly, the 
law of decrease of the electron temperature became universal at large times, 
i.e. did not depend on the initial conditions. (For example, when the initial 
temperatures T e varied by 30 times, their difference after a few microseconds 
was only 2-i-3 times and decreased further in the course of the subsequent 
evolution |8J.) Secondly, which is even more interesting, the asymptotics was 
measured to be T e oc f-C 1 - 2 ^ - 1 ) « t^ 1 |9J instead of t~ 2 , which would be ex- 
pected for an ideal rarefied gas without internal degrees of freedom (7 = 5/3) 
at the inertial stage of its expansion (i.e. when the plasma cloud expands with 
a constant velocity, so that its characteristic size increases linearly with time, 
R oc t). 

The most straightforward way to interpret the substantially slower decrease in 
the electron temperature is to take into account a heat release by the recom- 
bination of charged particles. In the particular case of atomic ions, the most 
efficient channel should be three-body recombination A + + e + e — > A + e, when 
one electron is captured by the ion, while the second electron carries away the 
excessive energy. Unfortunately, the recent attempts of numerical simulation of 
the observed law of temperature evolution with the three-body recombination 
did not lead to the satisfactory results; see, for example, Fig. 3a [5]) 

The aim of our paper is to show that all the experimentally observed features 
of temperature evolution (namely, both the establishment of the universal 
asymptotics, independent of the initial conditions, and its particular form, 
close to t^ 1 ) can be naturally explained by the model of virialization of the 
charged particle velocities, i.e. actually due to changing the equation of state 
of the ultracold plasma under the presence of strong electron-ion correlations. 

1 Let us mention that the method of plotting the temperature curves in Fig. 3a [9] 
is slightly confusing. According to the physical sense of the problem, the various 
laws of evolution should be compared at the same initial temperature; while in the 
above-mentioned figure they were presented at various initial temperatures. As a 
result, at first sight, the disagreement between the curves appears at small rather 
than at large times. 
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Thereby, in the first approximation, it is not necessary to take into account 
any inelastic processes, such as heat release by the three-body recombination. 



2 Theoretical Model 



Our analysis will consist of the three basic steps. First of all, let us consider 
a sufficiently small (but macroscopic) element of the expanding plasma where 
thermodynamic equilibrium is supposed to be established. Then, we can de- 
scribe it by the multi-particle distribution function of the following general 
form: 



/(r e i,...,r eJVe ,v el ,...,v e7Ve ) = A/expj-^-^ 
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where r en and v en are the coordinates and velocities of electrons, r in are the 
coordinates of ions, and Af is the normalization factor. (Kinetic energy of the 
ions is ignored here, since it is usually much less than the kinetic energy of 
electrons in the particular experimental setups. In principle, the same consid- 
eration can be conducted with the kinetic energy of ions if the thermodynamic 
equilibrium is established between all kinds of the charged particles.) 

Despite a very complex form of the potential energy U in the regime of strong 
Coulomb's interaction, the average value of some quantity F depending only 
on the velocities v en can be calculated quite easily: 
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because the integrals / exp{— U(r e , r,i)/k B T e } dr e in the numerator and de- 
nominator exactly cancel each other. (Here, v e , r e , and designate the entire 
sets of velocities and coordinates of the electrons and ions, respectively.) 

Therefore, the average kinetic energy per one particle will be 



(k) = (3/2) k B T e . 



(3) 



This formula looks exactly as for the ideal gas, but it is actually applicable for 
a plasma with arbitrary strong Coulomb's interaction U between the particles. 
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Next, at the second step of our consideration, the average kinetic energy (k) 
can be related to the average potential energy (u) by the virial theorem for 
the Coulomb's field [TO] S 



(k) = (1/2) \(u)\, 



(4) 



which is also valid at the arbitrary strength of interparticle interaction. Strictly 
speaking, the virial theorem is applicable only to the system of particles ex- 
periencing a finite (i.e. restricted in space) motion. Nevertheless, we can ex- 
pect that this theorem will be sufficiently accurate also for a freely expanding 
plasma cloud if the characteristic time of variation in its macroscopic param- 
eters (> 1CT 5 s for the particular experimental setup PTTT] ) is considerably 
greater than the microscopic periods of motion of the electrons (10~ 9 -^lCr 7 s). 

Finally, at the last step of our consideration, the average potential energy 
can be evidently expressed through the characteristic distance between the 
particles: 



Therefore, combining all the above formulas, we get T e oc n 1 ^ 3 . In particular, 
if the cloud expands inertially (i.e. linearly in time) and, consequently, its 
concentration changes as t~ 3 , we get: 



In summary, we conclude that the presented model of virialization of the 
charged particle velocities well explains both experimental features of ther- 
mal evolution of the plasma, namely: (a) the system "forget" in the course of 
time about its initial temperature, i.e. the plasma clouds with various initial 
temperatures evolve by the same way; and (b) the particular functional de- 
pendence is close to t _1 instead of the intuitively expected t~ 2 . The inelastic 
processes, such as three-body recombination, are not of importance in this 
model. 



3 Numerical Simulations 

To verify the above theoretical estimates, we performed ah initio simulation 
of the plasma dynamics, based on the numerical solution of the equations of 



Since the virial theorem directly relates the energies averaged over time, we need 
to assume also that the system is ergodic, i. e. the quantities averaged over time are 
equal to the ones averaged over ensemble. 




(5) 
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classical mechanics for the multi-particle system. Our approach differed from 
the earlier works in the following aspects. 

First of all, the authors of most of the previous simulations of ultracold plasmas 
tried to include into consideration as many particles as possible. Unfortunately, 
to reduce the computational cost, they had to use some simplifications of the 
equations of motion of the light particles (electrons), such as the particles-in- 
cell (PIC) method [12] or Vlasov approximation for the electrons [13]. Both 
these approaches are based on the introduction of average electric fields and, 
therefore, completely ignore strong individual electron-ion interactions (large- 
angular scattering), which are just responsible for the virialization of velocities. 
As distinct from these approximations, we did not intend to simulate as many 
particles as possible but tried to integrate the equations of motion of the 
electrons in the real electric microfields as accurately as possible, without 
using any extra simplifications. (The ions were assumed to be very heavy and 
moved by a purely inertial law.) 

Yet another well-known problem in modeling of the expanding plasmas is a 
considerable change of the spatial scale of the system (and, therefore, the 
amplitude of Coulomb's forces) during its evolution. As a result, it is quite 
difficult to choose the method of numerical integration ensuring a reasonable 
accuracy in the entire time interval. We resolved this problem by introducing 
a "scalable" coordinate frame, expanding in space with the average velocity 
of plasma outflow. In other words, the initial equations of the electron motion 

d 2 v ek /dt 2 = F ek , k = 1, . . . , N e (7) 

(where F ek is the total Coulomb's force acting on /c'th electron from all 
other electrons and ions) after the introduction of dimensionless variables 
r* = r/l, t* = tjr and transformation to the coordinate frame expanding 
with plasma, I = l (1 + u^t*), can be reduced to 

r ek +2u* (i+u* 1*)- 1 Kk = (i + «S 0~ 8 K k ■ (8) 

Here, r = (m/Ze 2 ) 1 ^ 2 Iq 3 ^ is the characteristic plasma time (on the order of 
the inverse Langmuir frequency), I is the characteristic distance between the 
particles (Iq is its value at the initial instant of time); Uq is the dimension- 
less velocity of the inertial expansion of the plasma cloud, determined by the 
relation u* = Uqt/L , where L is the size of the computational cell used in 
our simulations, and u is the velocity of its boundary; F* fc are the dimen- 
sionless Coulomb's forces, and the dot denotes derivative with respect to the 
dimensionless time t*. 

Therefore, as follows from equations ([8]), the effect of inertial plasma outflow 
in the expanding coordinate system looks like an effective dissipative force, 
proportional to the electron velocities. Consequently, the temperature of the 
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Fig. 1. Simulated electron temperature as function of time in the logarithmic scale. 
The inclined dashed lines show the power- like dependences T e oc t a \ the values of a 
being presented near the respective lines. 



electron gas is determined by the competition between two effects: (a) acceler- 
ation and heating of the electrons due to Coulomb's interactions with ions and 
(b) their deceleration and cooling by the above-mentioned dissipative forces. 



It is important to emphasize also that transformation from the equations (j7j) 
to (jSj) enabled us to perform a numerical integration in the fixed region of 
dimensionless coordinates and, therefore, to avoid the problem of quick losing 
the computational accuracy when the spatial scale of the system changes very 
much. 



The results of our simulations are presented by crosses in Fig. HJ If we do 
not take into account the earliest time interval, when the relaxation processes 
occur, the computational points in logarithmic coordinates are located almost 
along a straight line, corresponding to the power law of evolution of the elec- 
tron temperature: T e ~ t a . The corresponding exponent was estimated to be 
within the limits a = —(1.08 -j- 1.25). This is quite close to the value a = — 1, 
following from the simple virial estimate, and is in excellent agreement with 
the experimental values a = —(1.1 -r- 1.3) [9]. 
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4 Conclusions 



Analysis of the available data on the evolution of the electron temperature in 
the freely expanding ultracold plasma clouds shows that a heat release due 
to three-body recombination cannot explain quantitatively the experimental 
results. On the other hand, taking into account the strong electron-ion cor- 
relations and the resulting virialization of the charged particle velocities (i.e. 
modification of the equation of state of the plasma) leads to the perfect agree- 
ment with the experimental data. 
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